The Quantum Separability Problem for Gaussian States 
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Determining whether a quantum state is separable or entangled is a problem of fundamental 
importance in quantum information science. This is a brief review in which we consider the problem 
for states in infinite dimensional Hilbert spaces. We show how the problem becomes tractable for a 
class of Gaussian states. 
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I. INTRODUCTION 



The concept of entanglement arose with the question of completeness of quantum theory . Nowadays entanglement 
I""!' is regarded as a fundamental property of certain quantum states and it appears to be an important physical resource. 
In some sense, entanglement is synonymous of inseparability because entangled states possess some global properties 
. that cannot be explained in terms of only the parties (subsystems) of the system. Roughly speaking, entangled states 
^ ' possess "strong" correlations among parties that cannot be explained within any classical local theory (because these 
O ■ would imply an instantaneous action at distance). Separable states may also exhibit correlations among parties, but 

these are purely classical and local, hence "weaker" than those underlying entanglement, 
^vq , Recently, the role of entanglement became important and often necessary in many different contexts like quantum 
^ ■ algorithms, quantum communication protocols, quantum cryptography, etc. (see e.g. 0)- So, the problem of deciding 
' whether a given quantum state is separable or entangled has become of uppermost importance. This can be called the 
> Quantum Separability Problem (QSP). Essentially, it represents an instance of a combinatorial optimization problems 
' called the Weak Membership Problem 0]. 

' Although there exists a number of characterizations of separability, there is still no feasible procedure to solve QSP 
in its generality (see e.g. Q and references therein). Concerning its computational complexity, QSP is a "difficult" 
problem. In fact, QSP has been proved to be NP-hard @. However, if we restrict ourselves to specific classes of 
quantum states, there are examples in which QSP can be efficiently solved. For instance, this is the case of states in 
Hilbert space of dimension 2 or 3 Q and certain finite sets of states Q. 

In infinite dimensional Hilbert spaces, Gaussian states give rise to an important class of states for which QSP 
is "easy" (see e.g. 013 and the reference therein). In this paper, we review the formulation of QSP for infinite 
dimensional Hilbert spaces and we show how to tackle the problem for the class of Gaussian states. 
$^ ' The paper is organized as follow. In Sec^ we review some basic notions of Quantum Theory. In Sec. IIIII we 
5^ , formalize the QSP. In Sec II VI we introduce the Gaussian states. In Sec|V]we develop a criterion for separability of 
Gaussian states. Finally, some conclusions are drawn in Sec lVIl 



II. BASIC NOTIONS OF QUANTUM THEORY 



In this section, we introduce some terms and notions of Quantum Mechanics needed to approach the paper. Of 
course, the expert reader may skip this section. 

In its standard formulation. Quantum Theory takes place in Hilbert spaces A Hilbert space 7i is a vector 

space over the field of complex numbers C endowed with an inner product (which induces a norm), that can have 
finite or infinite dimension. We use the so-called Dirac notation for a vector 1-0). Its dual is (-01. Then, the inner 
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product between two states \'4)) and \4>) reads {'4'\(j)) € C. The norm of a vector |^) results \\ \'4))\\ = ^JJ^p\^pj. The 
following two postulate fix the mathematical representation of quantum states: 

Postulate II. 1. The space of states of a physical system is a Hilbert space. The states are described by unit norm 
vectors in such Hilbert space. 

Postulate II. 2. The space of states of a composite system is the tensor product of Hilbert spaces of subsystems. 

The structure of Hilbert space naturally leads, when considering composite systems, to the concept of entanglement. 
In fact, there exist states of the whole system that cannot be factorizcd into states of the subsystems. 

Example II. 1. Let \'ip)i^ \ipi_)i be two orthogonal states in Tii and \rfj2, |?7±)2 ^'^^ orthogonal states in Ti.2. Then, 
® \v)2 ^ '^1 ^ "^^2 as well as {a\ilj)^ (g) \r])2 + b "8) \V-i-)2) G '^1 ® "^^2, with a, 5 £ C. The first can be factorizcd 

into states of the subsystems; this is not the case for the second one. 

It is fashinating that this seemingly abstract mathematical notion has a large impact in the description of the 
quantum mechanical world. 

The above postulates can be generalized in terms of mixture of states, {pj, \tpj)}, where denotes for the probability 
for the system to be in the state {tpj). This can be done by introducing the notion of density operator: 

Definition II. 1. A density operator p is a non-negative, self-adjoint, trace-one class operator which is also positive 
semi-definite (that is {■ilj\p\^j) > y\'ip) G Ti.. 

Thus we can represent the mixture {pj, \ipj)} by the density operator p ~ J2jPjH'j ){''Pj\- 

Definition II. 2. A state p of a composite bipartite system is said to be separable iff it can be written in the form 

with non-negative pj 's such that '^jPj = 1; o-nd where pf \ pf^ are density operators of the subsystems; the state is 
said to be entangled otherwise. 

The physical quantities of a system that can (in principle) be measured are called observables. The next postulate 
fixes the mathematical representation of observables: 

Postulate II. 3. To physical observables correspond self-adjoint operators. The possible measurement results on the 
observable O are the eigenvalues of the associated self-adjoint operator O. The expectation value is (O) = Tr(Op). 

Restrictions on expectation values are imposed by the following famous principle: 

Principle II. 1 (The Uncertainty Principle). Any two observables A and B in H must satisfy, for all quantum 
states, the following inequality: 

((AA)2)((AB)2)>i|([A,B])|', (2) 
where AO = O — (O) and [A, B] = AB — BA is the commutator. 



III. THE QUANTUM SEPARABILITY PROBLEM 



In this section we introduce the Quantum Separability Problem. Let us consider a quantum system with two parties 
associated to a Hilbert space Hi (81 7^2 — (EiC^. Notice that such a Hilbert space is isomorphic to M*^ ^ and it is 

endowed with the Euclidean inner product {X , Y) = Ti{XY) which induces the corresponding norm = '\J tr(X^) 

and distance measure \\X — Y\\. Let V C Hi (8) H2 denote the set of all density operators. The set of bipartite 
separable quantum states, S C T>, is defined as the convex hull of the separable pure states 8) |''7)2(f?|} where 

\ip}i (resp. \ri)2) is a normalized vector in C*^ (resp. C^). An arbitrary density matrix in I? is parametrized by 
M'^N'^ — 1 real variables. Since we deal with continuous quantities, in defining the separability problem wc cannot 
allow infinite precision, so we need to introduce a precision parameter S G M+. 

Definition III.l (The Quantum Separability Problem). Given p ^ T> and a precision 5 assert either p is: 
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• Separable; there exists a separable state a such that \\p — ^\\ < ^ 
or 

• Entangled: there exists an entangled state f such that j|/5 — fH < |. 

In this formulation, this problem is equivalent to an instance of a combinatorial optimization problem called Weak 
Membership Problem In its complete generality, QSP has been shown to be NP-hard |^. Thus, any devised 
test for separability is likely to require a number of computational steps that increases very quickly with M and 
N. For MN < 6 the positivity under Partial Transpose (see the next section) represents a necessary and sufficient 
test 0. Otherwise, there only exist sufficient 'one-sided' tests for separability. In these tests, the output of some 
polynomial-time computable function of p can indicate that this is certainly entangled or certainly separable, but not 
both (see e.g. Q and reference therein). 



IV. GAUSSIAN STATES 



In this section we introduce Gaussian states. Let us now move to il/, N s- oo, thus considering two infinite 
dimensional Hilbcrt spaces Tii and 7^2- In such spaces we can introduce continuous spectrum self-adjoint operators 
corresponding to canonical position and momentum variables llO^ . Let us arrange them into four-dimensional column 
vectors 

v"^ = (91, Pi, 92, P2), = {xi,yi,X2,y2)- 
The operators in v obey commutation relations [To| that take the compact form 

[vq, V/3] = a,/? = 1,2,3,4, (3) 

with 

There is a one-to-one correspondence between density operators and c-number Wigner distribution functions in 
phase space the space of variables z, i.e. R*, in this case. 

Definition IV. 1. For a given density operator p in Hi and H2 the corresponding Wigner function is defined as 
followsmj 



W{z) Tr (pf (z)) , (5) 



where 



In turn, it results 



1 



T{^) I d'^' cxp zz'^ • (v - z) . (6) 



p = j d**zl^(z)r(z). (7) 

A density operator p has finite second order moments if Tr{pq^) < 00 and Tr(/5p|) < 00 for all j. In this case we 
can define the vector mean m as 



m :=Tr(/5v) 

= J d^zzW{z), (8) 



and the real symmetric correlation matrix V as 



i({Av„,Av0}), a,/? = 1,2, 3, 4, (9) 
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where {Avq,, Av/3} = Avq,Av^ + Av/jAv^ is the anticommutator. It results 



Tr(p-{Av„,Av^} 



j d'^z (z - m)„ (z - m)^ W{z) 



A given V is the correlation matrix a physical state iff it satisfies 



K = V + -n > 0, 



(10) 



(11) 



as consequence of the Uncertainty Principle [3 and commutation relation (jSJ . The correlation matrix forms a 4 x 4 
matrix that transforms as an irreducible second rank tensor under the linear canonical (symplectic) transformations 
and has 4 invariants. If we write the correlation matrix in the block form 



V = 



A C 
B 



(12) 



the invariants are detA, dciB, detC and Tv{AJCJBJC'^J). The condition ^ implies A > 1/4 and B > 1/4. 
Moreover, Ea. (|ll|l can be read as 



det AdetB- Tr{AJCJBJC^J) - - (det A + det B) 



- - det C) > 0. 
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It is also worth remarking that any correlation matrix can be brought into the standard form 



V 



(a c \ 

a d 
c b 



V 



(13) 



(14) 



with a,b,c,d G K, by effecting suitable local canonical transformations corresponding to some element of Sp(2,R)x 
Sp(2,R) C Sp(4,R). Now we are ready to give the definition of Gaussian state: 

Definition IV. 2. A state p is called Gaussian if its Wigner Junction takes the form 



W{7.) 



1 



: exp 



-i(z-m)^V""^(z-m) 



(15) 
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with m a real A-vector and V a real symmetric 4x4 -matrix. 

One can show that m is indeed the mean and V is the correlation matrix. These define the Gaussian state uniquely. 
In what follows, we simply consider the case m = 0, because m can be easily removed by some local displacement 
and thus has no influence on the separability or inseparability of the state. 

V. A SEPARABILITY CRITERION FOR GAUSSIAN STATES 



In this section, we describe how to solve QSP for Gaussian states. Let us consider a separable state psep of the 
form CQl in the Hilbert space T^i (8) • Let us choose a generic couple of observables for each subsystem, say fj , Sj 
on Hj (j = 1,2), with 



Cj^i[fj,Sj], j = l,2. 
Then, we introduce the following observables on Hi <E) 7^2 '■ 

u = aifi + a2f2 , 

V = biSl + b2S2 , 

with Qj, bj e M. From the the Uncertainty Principle |21 it follows that every state p on Hi ^ H2 must satisfy 



{iAur){{Avr)> 



\aibi{Ci) + 0262 (C2) 



(16) 



(17) 



(18) 



However, for separable states, a stronger bound exists. We have in fact the following theorem [iz 
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Theorem V.l. For any separable state the following implication holds: 

Psep =^ {{^uf){{Mf)>W\ 

where 

1 



with 



W=-(|ai6i|Wi + 102621 W2) 



W.EE^Pfc |(C;-),|, J = 1,2, 



(19) 



(20) 



(21) 



being {C,)k = 'Tv[C,jj^\ 
The theorem can be proved with the help of a family of linear inequalities 

a((Au)2) +/3((Aw)2) > 2v/^W, 



(22) 



which must be always satisfied by separable states. The convolution of such relations gives the condition l|19|l. 
representablc by a region in the ((Ati)^), ((A-0)^) plane delimited by an hyperbola. 
Notice that, since 



(23) 



the following inequalities hold 



W>-(|ai&i| |(Ci)| + |a2&2| |(C2)| 
>i(|ai&i (Ci)+a2&2 (C2)|). 



(24) 



In particular, Eq. H24|l tells us that the bound (|19|l for separable states is much stronger than Eq. H18|) for generic 
states. Moreover, Eq. H24() gives us a simple separability criterion. In fact, while W is not easy to evaluate directly, 
as it depends on the type of convex decomposition that one is considering, the right hand side of Eq. H24|) is 
easily measurable, as it depends on the expectation value of the observables Cj. Then, we can claim that Ea. (|19|l is 
a necessary criterion for separability, i.e. 



{{^uf){{Mf) < 



p entangled . 



(25) 



Example V.l. An important simplification applies when the observable Cj is proportional to the identity operator, 
e.g. fj = Qj and Sj = pj. In such a case, Ea. ^l8\) reduces to 



1 



((A«)^)((At>)^)>-, 



(26) 



while Eg. i^iy^ reduces to 

((Au)2)((Ai))2)>l, (27) 
Let us now consider the case in which fj, Sj are linear combinations of canonical observables qj and pj, i.e. 

0.3 



ri = qi 



ai 



-Pi 



04 

r2 = 92 H P2 

02 



si = Pi + —qi 
bi 



S2=P2 + -r-l2 , 

02 



(28) 



where 03, a4, 63, 64 G K are generic real parameters. Then Eq. 119() . taking into account Eq.Q, becomes 



{{Au)^){(Avf) > -( laifei - 0363! + |a262 - 0464! 



(29) 
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that should be compared with 

{{Au)^) + ({Avf) > \aibi ~ asbsl + \a2b2 - aM . (30) 

It is easy to verify that, given aj.bj [j ~ 1,. . . ,4), the "product condition" (|29|l impUes the "sum condition" (|3U|I . 
However, if we require Eqs. (|29() and H30() to be verified for all possible values of the coefficients aj,bj, the two are 
equivalent since it is possible to re-obtain one from another using a convolution trick, like the used with Eas. H19|l and 
(|22|l (the one-to-one correspondence between quadratic and linear tests under all circumstances has been also pointed 
out in Ref.O). 

It turns out that the restriction 

{{Auf) + {{Avf)>\aibi-a3b3\ + \a2b2-aibil Vaj-,6jeR, (31) 

is necessary and sufficient for separability of Gaussian states [T^ ITsI l . 

However, solving QSP by testing the condition (|31|l would be hard from a complexity point of view, due to the 
presence of the universal quantifier at right hand side. 

Nevertheless, the condition H31|) can be rephrased in a simpler way. First notice that the uncertainty relation 
satisfied by all (separable and inseparable) states 

((A?/)2)((A«)2) > i|ai6i - 0363 + a2&2 - (32) 

and corresponding to it and v formed from Ea. p9() . is equivalent to 

((Am)2) + ((Au)2) > |ai5i - 0363 + 0262 - aibil, (33) 

as much as like Eas. (|^ and (ISOJ. 

Then, what is the relation between conditions H3()(l and H33(l ? 
They are simply related by the partial transpose transform 

PT: V — >Av, A = diag(l, 1,1,-1). (34) 

This operation inverts P2, leaving qi, pi, and 92 unchanged [Toj. In fact, separable states satisfies the usual uncertainty 
relation H33() and the analogous one obtained under partial transpose; thus these satisfy the condition (|30|1 . 

On the other hand, the transformation l|34|l changes the correlation matrix asV^V = AT^A. Hence, the compact 
uncertainty relation Hll|l becomes 

V+^Q>0. (35) 

Expressed in terms of invariants, the condition (|35|l for V takes a form identical to (|13|l . The signature in front of 
det C in the second term on the left hand side is changed. Thus, if we write 

/I 

fiV) :=detAdetS+ f - - |detC| 

-ti{AJCJBJC'^J) - ^{dctA + detB), (36) 

the requirement that the correlation matrix of a separable state has to obey ()35|l . in addition to the fundamental 
uncertainty principle (|11|) . can be stated as follow 

Theorem V.2. A bipartite Gaussian state is separable iff f{V) > 0. 

The necessity follows from theorem lV.ll The sufficiency follows from the fact that Gaussian states with correlation 
matrix having det C > are separable [lj| . 

The statement IV. 21 is equivalent to the condition (p?T|l . but much more effective to be used. 

Given the standard form H14|l of the correlation matrix V we can consider the space of all possible Gaussian states 
as isomorphic to W^, while the set of physical states is a subspace ^ C IR^ defined through Hll|l . Furthermore, the 
equation f{V) — reads 

f{V) = 4(a6 - c^){ab ~ (f) - {a^ + b^) - 2\cd\ " ^ = 0. (37) 
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The equation defines the surface S of the subset S C G oi separable states. Then, by simply evaluating / we can 
say whether a given state (point in Q) is within S (hence separable) or not (hence entangled). This is an easy 
computational task that can be efficiently accomplished. In reality, taking into account a finite accuracy (5, wc can 
only say that the state is almost separable (resp. almost entangled) within S. Nevertheless, if we want to assert that 
the state is strictly separable (resp. strictly entangled), we have to be sure that the distance of the state p from the 
surface S is greater than 1/6. That is 

min||/}-/}'|| > i. (38) 
pes 

According to Scc lIIII the distance between two states is considered as ||p — = ^Tr [{p — p')"^] and for Gaussian 
states this can be expressed through Wigner functions (hence correlation matrices) as 

\\p^p'\\= ld^z[W{z)-W'{z)f. (39) 

Such a task can be efficiently accomplished with the aid of geometrical arguments and simple algorithms. For instance, 
a software package that efficiently find all hyperplanes tangent to the surface S, from which evaluate the l.h.s. of 
Ea H38() , is already available . 



VI. CONCLUSIONS 



Summarizing, we have given a brief review of QSP for Gaussian states of two parties. The problem has been 
approached by developing tests that involve variances to arrive at an efficient solution based on the invariance (pos- 
itivity) of only separable states under partial transpose. Notice that this argument can be further generalized to 
partial scaling transforms to which partial transpose belongs. In fact, while K and V are always invariant under lin- 
ear canonical transformations, they are not invariant under scale changes on the v that are not contained in Sp(4, R). 
In particular under partial scaling K is not necessarily positive definite [T^ . These arguments could be extended to 
multipartite systems, with e.g. N degrees of freedom. Starting from the uncertainty relation K = V + ^il > 0, we 
can perform an arbitrary scaling described by the real vector x = {xi, X2, ... X2n) and then compute 

^.x = yx _^ * yx ^ AxVAx, (40) 

with Ax = diag(a;i, a;2, ...X2n)- The 2N real quantities x parameterize the Abelian scaling semigroup with the 
requirement that 

|a;iX2| > 1, 1x3X4! > 1, . . . , |a;2jv-iX2Ar| > 1. (41) 

The necessary condition for the separability of the state is 

> 0, Vx. (42) 

Notice, however, that for multipartite systems besides separability (resp. inseparability) there can be the possibility 
of partial separability (resp. partial inseparability), e.g. separability of a subsystem with respect to the others which 
in turns are entangled Hence QSP becomes much more subtle and even for Gaussian states it is not completely 
understood. 
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